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1. Introduction and Overview 

These lectures focus on three topics. The first is an overview of the current paradigm of 
early universe cosmology, the inflationary universe scenario. The second topic is a pedagogical 
presentation of the theory of cosmological perturbations, the main tool of modern cosmology which 
allows us to connect theories of the very early universe with observational data. The third topic is 
string gas cosmology, an attempt to construct a new scenario of the very early universe based on 
established principles of string theory. 

Over the past two and a half decades, cosmology has become a science dominated by data of 
rapidly increasing accuracy. Today, we have three-dimensional maps of the distribution of galaxies 
in space which contain more than one hundred thousand galaxies [[I], |2|]. They clearly indicate that 
luminous matter in the universe is neither uniformly nor randomly distributed. There are clear 
patterns to be seen: clusters of galaxies, superclusters, filaments and voids (regions of space empty 
of galaxies). The distribution can be quantified in terms of the luminosity power spectrum. A key 
challenge for cosmology is to understand the origin of these patterns in the distribution of matter. 

Another observational window in cosmology is the cosmic microwave background (CMB) 
radiation. Overall, this radiation is characterized by a surprising isotropy. At a fractional level of a 
bit less than 10~ 4 , however, there are anisotropics. These can be quantified in terms of their angular 
power spectrum. First, the sky map (two-dimensional) of anisotropics is expanded in spherical 
harmonics Y\ m : 

Ar °° l - 

^r(0,<P) = E £ a lm Y lm (e,(p) (1.1) 

1 l=l m= -l 

where 6 and (j) are the usual angles on the surface of the sphere. If the fluctuations are due to a 
Gaussian random process with no distinguished direction in the sky, then the complete information 
about the fluctuations is given by the ensemble average (denoted by pointed parentheses) of the 
coefficients a/ m : 

ci =< \ai, n \ 2 >, (1.2) 

These c/ coefficients define the angular power spectrum of CMB anisotropics. Figure 1 is the 
full sky map of CMB anisotropics from the WMAP experiment [||]. Figure 2 shows the resulting 
angular power spectrum )(what is plotted on the vertical axis is /(/ + l)C//(27r)). The key features 
are the flat region at small values of I (large angular scales) and the characteristic oscillations of 
the spectrum at intermediate angular scales. Another crucial challenge for cosmology is to explain 
both the overall isotropy of the CMB, and the specific patterns of anisotropics. 

According to our present understanding, we must look to the very early universe to find an 
explanation for the observed structures. The reason is that in Standard Big Bang cosmology (SBB), 
which well describes the cosmological evolution at late times (times later and including the period 
of nucleosynthesis) the physical wavelength of fixed comoving scales is increasing less fast than 
the Hubble radius (an important length scale which is defined and whose role is described at the 
end of this section). Hence, the scales which are currently observed were outside of the Hubble 
radius at early times, and no causal structure formation scenario is possible . 

'Topological defect models |Q, ^, []] provide a way to circumvent this reasoning. They, however, also involve new 
physics of the very early universe. 
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Figure 1: All sky map of the temperature anisotropies in the CMB from the WMAP satellite experiment. 
Credit: NASAAVMAP Science Team. 

Angular Scale 




Multipole moment (I) 

Figure 2: The power spectrum of CMB anisotropies as computed from the WMAP satellite experiment. 
The horizontal axis is the angular quantum number I, the vertical axis gives the power of the CMB on the 
respective scales. The dots represent the data points (with their errors bars indicated), the solid curve is the 
prediction of the best-fit inflationary model. Credit: NASAAVMAP Science Team. 



It is a remarkable success of inflationary cosmology, our current paradigm of early universe 
cosmology, that it, in addition to explaining why the universe is large, spatially flat, and containing 
a large amount of entropy, provides a causal mechanism for the origin of inhomogeneities in the 
universe. The solid curve in Figure 2 represents the predictions of inflationary cosmology 2 The 

2 Several cosmological parameters, e.g. the current value of the cosmological constant and the fractional baryon 
density, have to be fixed in order to obtain this excellent agreement. The number of free parameters, however, is much, 
much smaller than the number of data points. 
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predictions were made more than 15 years before the data shown in Figure 1 were collected. 

In these lectures, however, I would like to focus less on the phenomenological successes of 
the inflationary paradigm, but more on the conceptual problems which current realizations of in- 
flationary cosmology are confronted with. These problems call out for the development of a new 
paradigm, a paradigm which must include new fundamental physics, a theory which better de- 
scribes space, time and matter and the highest densities. The best current candidate for such a 
theory is superstring theory. Hence, in the third main part of these lectures (Section 4) I will ex- 
plore the possibility that string theory may lead to a new paradigm of early universe cosmology. 
My approach to string cosmology in this section is complementary to the one taken in the lectures 
of Cliff Burgess at this school (see [^]), in which avenues of obtaining inflation in the context of 
models coming from string theory are explored (see also []8j, ^] for other reviews of such avenues). 

The theory of cosmological perturbations plays a key role in modern cosmology, since it pro- 
vides the techniques with which to calculate, in the context of any given scenario of the very early 
universe, the generation and evolution of the predicted inhomogeneities in the matter distribution 
and anisotropics in the CMB, and thus allows for a comparison between fundamental theory and 
observational data. The second part of these lectures (Section 3) provides a pedagogical overview 
of this theory. 

I begin, however, with a discussion of inflationary cosmology, the current paradigm of cos- 
mology. 

To establish our notation and framework, we will be taking the background space-time to be 
homogeneous and isotropic, with a metric given by 

ds 2 = dt 2 -a(t) 2 dx 2 , (1.3) 

where t is physical time, dx 2 is the Euclidean metric of the spatial hypersurfaces (here taken for 
simplicity to be spatially flat), and a(t) is the scale factor. The scale factor determines the Hubble 
expansion rate via 

H(t) = ~(t). (1.4) 

The coordinates x used above are comoving coordinates, coordinates painted onto the expanding 
spatial hypersurfaces. 

We are interested in tracking the time evolution of the physical wavelength of the currently 
observed patterns in the distribution of matter and of the CMB anisotropies. Since the patterns are 
assumed to be frozen in in comoving coordinates, the physical wavelength scales as a(t). 

A key length scale in cosmology is the Hubble radius 

l H (t)=H-\t), (1.5) 

defined to be the inverse Hubble expansion rate. As will be explained later, the Hubble radius is the 
maximal distance that microphysics can act coherently over a Hubble expansion time - in particular 
it is the maximal distance on which any causal process could create fluctuations. 

2. Inflationary Cosmology: The Current Paradigm 

Standard Big Bang (SBB) cosmology, the precursor to inflationary cosmology as the paradigm 
for the evolution of the early universe, is based on a classical physics description of both space-time 
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Figure 3: Space-time diagram (sketch) showing the evolution of scales in inflationary cosmology. The 
vertical axis is time, and the period of inflation lasts between f; and tg, and is followed by the radiation- 
dominated phase of standard big bang cosmology. During exponential inflation, the Hubble radius H^ 1 is 
constant in physical spatial coordinates (the horizontal axis), whereas it increases linearly in time after tR. 
The physical length corresponding to a fixed comoving length scale labelled by its wavenumber k increases 
exponentially during inflation but increases less fast than the Hubble radius (namely as J 1 / 2 ), after inflation. 



(via Einstein's theory of General Relativity) and matter (a superposition of two perfect fluids, the 
first describing pressureless matter - cold matter - the second describing radiation - the CMB). The 
key phenomenological success of the SBB model is the prediction of the existence and black body 
nature of the CMB, the black body nature of which was confirmed with spectacular accuracy by 



the COBE satellite [jig] and UBC rocket experiments Qll|]. 

On the other hand, the SBB scenario leaves many crucial questions un-answered Jl2|]. Why is 
the universe so close to spatially flat? Why is it so large and contains such a large entropy? Why is 
the CMB isotropic to an accuracy of better than 10~ 4 (after subtracting the dipole contributions due 
to our motion relative to the rest frame of the CMB and the effects due to the emission of our own 
galaxy)? Most importantly, what is the origin of the observed inhomogeneities in the distribution 
of matter and of the small CMB anisotropics? These are the "flatness", "entropy", "horizon" and 
"structure formation" problems of SBB cosmology. SBB also suffers from conceptual problems: 
the initial cosmological singularity tells us both that the theory must be incomplete, and that it is 
based on using the wrong fundamental physics input close to the singularity 3 . 



Inflationary cosmology [12] (see also [13, 14, 15]) provides an solution of the horizon, flatness 
and entropy problems. In addition, it provides a mechanism for the origin of structure in the 



universe based on causal physics Q16| ] (see also y 

The idea of inflationary cosmology is to assume that there was a period in the very early 
Universe during which the scale factor was accelerating, i.e. a > 0. This implies that the Hubble 
radius was shrinking in comoving coordinates, or, equivalently, that fixed comoving scales were 
"exiting" the Hubble radius. In the simplest models of inflation, the scale factor increases nearly 
exponentially. 



As we will see later on, these conceptual problems persist in inflationary cosmology. 
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As illustrated in Figure [3], the basic geometry of inflationary cosmology provides a solution 
of the fluctuation problem. As long as the phase of inflation is sufficiently long, all length scales 
within our present Hubble radius today originate at the beginning of inflation with a wavelength 
smaller than the Hubble radius at that time. Thus, it is possible to create perturbations locally using 
physics obeying the laws of special relativity (in particular causality). As will be discussed later, it 
is quantum vacuum fluctuations of matter fields and their associated curvature perturbations which 
are predicted to be responsible for the structure we observe today. 

Postulating a phase of inflation in the very early universe solves the horizon problem of the 
SBB, namely it explains how the causal horizon at the time t rec when photons last scatter can be 
larger than the radius of the past light cone at t rec , the part of the last scattering surface which is 
visible today in CMB experiments. Inflation explains the near flatness of the universe: in a deceler- 
ating universe spatial flatness is an unstable fixed point of the dynamics, whereas in an accelerating 
universe it becomes an attractor. Another important feature of inflation is that the volume of space 
increases exponentially at constant energy density. If this energy density is successfully converted 
to ordinary matter at the end of the period of inflation, then the entropy of the universe is expo- 
nentially larger after compared to before inflation. In addition, with the exponential expansion of 
space it is easy to produce a universe of our size today from a Planck scale universe at the initial 
Planck time, something which is not possible in the SBB model. Thus inflation explains the large 
current size and entropy of the universe. 

Let us now consider how it is possible to obtain a phase of cosmological inflation. We will 
assume that space-time is described using the equations of General Relativity 4 . In this case, 
the dynamics of the scale factor a(t) is determined by the Friedmann-Robertson- Walker (FRW) 
equations 

(^) 2 = SKGp (2.1) 

and 

^ = -47iG(p + 3p) (2.2) 

where for simplicity we have omitted the contributions of spatial curvature (since spatial curvature 
is diluted during inflation) and of the cosmological constant (since any small cosmological constant 
which might be present today has no effect in the early Universe since the associated energy density 
does not increase when going into the past). In the above, p and p denote the energy density and 
pressure, respectively. From (2.2) it is clear that in order to obtain an accelerating universe, matter 
with sufficiently negative pressure 

P < -~P (2-3) 

is required. Exponential inflation is obtained for p = — p. 

Conventional perfect fluids have positive semi-definite pressure and thus cannot yield inflation. 
However, we know that a description of matter in terms of classical perfect fluids must break down 
at early times. An improved description of matter will be given in terms of quantum fields. Scalar 
matter fields are special in that they allow at the level of a renormalizable action the presence 



4 Note, however, that the first model of exponential expansion of space Jl5| ] made use of a higher derivative gravita- 
tional action. 
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of a potential energy term. The energy density and pressure of a scalar field <p with canonically 
normalized action 5 

S = Jd*x y /=g[-d, l q>dl l q>-V(q>)] (2.4) 

(where Greek indices label space-time coordinates, V(<p) is the potential energy density, and g is 
the determinant of the metric) are given by 

P = \(<{>) 2 + l -a- 2 {V<p) 2 + V{<p) 

P = \(<P?-\a-\V<pf-V(<p). (2.5) 

Thus, it is possible to obtain an almost exponentially expanding universe provided the scalar field 
configuration 6 satisfies 

\(V P <?) 2 « V(q>), (2.6) 
l -{<p) 2 « V(<p). (2.7) 

In the above, V p = a~ l V is the gradient with respect to physical as opposed to comoving coordi- 
nates. Since spatial gradients redshift as the universe expands, the first condition will (for single 
scalar field models) always be satisfied if it is satisfied at the initial time 1 . It is the second condition 
which is harder to satisfy. In particular, this condition is in general not preserved in time even it is 
initially satisfied [20]. 

It is sufficient to obtain a period of cosmological inflation that the slow-roll conditions for q> 
are satisfied. Recall that the equation of motion for a homogeneous scalar field in a cosmological 
space-time is (as follows from (2.4)) is 

<p + 3H<p = -V'((p), (2.8) 

where a prime indicates the derivative with respect to (p. In order that the scalar field roll slowly, it 
is necessary that 

(f> < 3H<p (2.9) 



such that the first term in the scalar field equation of motion ( |2.8| ) is negligible. In this case, the 
condition (|L7|) becomes 

{—f < 487TG (2.10) 



and (2.9) becomes 

V" 

— < 24ttG. (2.11) 

In the initial model of inflation using scalar fields ("old inflation" [|l^]), it was assumed that 
cp was initially in a false vacuum with large potential energy. Hence, the conditions for inflation 



5 See jlil l for a discussion of fields with non-canonical kinetic terms. 
6 The scalar field yielding inflation is called the inflaton. 

7 In fact, careful studies [ |l^ ] show that since the gradients decrease even in a non-inflationary backgrounds, they can 
become subdominant even if they initially dominate. 
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are trivially satisfied. To end inflation, a quantum tunneling event from the false vacuum to the 



true vacuum y21[ ] was invoked (see e.g. [ ]22[ ] for a pedagogical review). This model, however, has 
a graceful exit problem since the tunneling leads to an initially microscopical bubble of the true 
vacuum which cannot grow to encompass our presently observed universe - the flatness problem 
of SBB cosmology in a new form. Hence, attention shifted to models in which the scalar field (p is 
slowly rolling during inflation. 

There are many models of scalar field-driven inflation. Many of them can be divided into three 



groups [23]: small-field inflation, large-field inflation and hybrid inflation. Small-field inflationary 
models are based on ideas from spontaneous symmetry breaking in particle physics. We take the 
scalar field to have a potential of the form 

V(q>) = l -X{<p 2 -o 2 )\ (2.12) 

where a can be interpreted as a symmetry breaking scale, and A is a dimensionless coupling con- 
stant. The hope of initial small-field models ("new inflation" [p4|]) was that the scalar field would 
begin rolling close to its symmetric point (p = 0, where thermal equilibrium initial conditions would 
localize it in the early universe. At sufficiently high temperatures, (p = is a stable ground state of 
the one-loop finite temperature effective potential W(<P) (see e.g. [ p2| ] for a review). Once the tem- 
perature drops to a value smaller than the critical temperature T c , (p = turns into an unstable local 
maximum of Vr(<p), and (p is free to roll towards a ground state of the zero temperature potential 
( [2.12 ). The direction of the initial rolling is triggered by quantum fluctuations. The reader can eas- 



ily check that for the potential ( |2.12 ) the slow -roll conditions cannot be satisfied if a <C m p i, where 



m p i is the Planck mass which is related to G. If the potential is modified to a Coleman-Weinberg 



[]25J] form 

A A r , I CO I K 1 



^) = ^>^4] +1 ^a 4 (2.13) 

(where a denotes the value of the minimum of the potential) then the slow-roll conditions can be 
satisfied. However, this corresponds to a severe fine-tuning of the shape of the potential. A further 



problem for most small-field models of inflation (see e.g. [gfj for a review) is that the slow-roll 
trajectory is not an attractor in phase space. In order to end up close to the slow-roll trajectory, 
the initial field velocity must be constrained to be very small. This initial condition problem of 
small-field models of inflation effects a number of recently proposed brane inflation scenarios, see 
e.g. [ p6| ] for a discussion. 

There is another reason for abandoning small-field inflation models: in order to obtain a suffi- 
ciently small amplitude of density fluctuations, the interaction coefficients of cp must be very small 
(this problem is discussed in detail at the beginning of Section 4). In particular, this makes it incon- 
sistent to assume that cp started out in thermal equilibrium. In the absence of thermal equilibrium, 
the phase space of initial conditions is much larger for large values of (p. 



This brings us to the discussion of large-field inflation models, initially proposed in [27] under 
the name "chaotic inflation". The simplest example is provided by a massive scalar field with 
potential 

V((p) = ^m 2 (p 2 , (2.14) 
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where m is the mass. It is assumed that the scalar field rolls towards the origin from large values 



of \<p\. It is a simple exercise for the reader to verify that the slow -roll conditions ( |2.10[ ) and ( J2. 1 1| ) 
are satisfied provided 

\<p\>^m pl . (2.15) 

Values of | <p | comparable or greater than m p i are also required in other realizations of large-field 
inflation. Hence, one may worry whether such a toy model can consistently be embedded in a re- 
alistic particle physics model, e.g. supergravity. In many such models V((p) receives supergravity- 
induced correction terms which destroy the flatness of the potential for |<p| > m p i. However, as 



discussed e.g. in [28], if the flatness of the potential is protected by some symmetry, then it can 
survive inclusion of the correction terms. As will be discussed later, a value of m ~ 10 13 GeV is 
required in order to obtain the observed amplitude of density fluctuations. Hence, the configuration 
space of field values with |<p| > m p i but V ((p) < m 4 pl is huge. It can also be verified that the slow -roll 
trajectory is a local attractor in field initial condition space [19], even including metric fluctuations 



at the perturbative level [ggj. 

With two scalar fields it is possible to construct a class of models which combine some of the 
nice features of large-field inflation (large phase space of initial conditions yielding inflation) and 
of small-field inflation (inflation taking place at sub-Planckian field values). These are models of 



hybrid inflation [30]. To give a prototypical example, consider two scalar fields q> and % with a 
potential 

V(q>,z) = \^x(X 2 -o 2 ) 2 + \™ 2 <P 2 -\g 2 <P 2 X 2 - (2-16) 

In the absence of thermal equilibrium, it is natural to assume that |<p| begins at large values, values 
for which the effective mass of % is positive and hence % begins at % = 0. The parameters in the 
potential ( |2.16i ) are now chosen such that <p is slowly rolling for values of |<p| somewhat smaller 
than m p i, but that the potential energy for these field values is dominated by the first term on the 



right-hand side of ( |2.16| ). The reader can easily verify that for this model it is no longer required 
to have values of | q>\ greater than m p \ in order to obtain slow -rolling 8 The field (p is slowly rolling 
whereas the potential energy is determined by the contribution from %. Once \ <p\ drops to the value 

|<p c | = (2.17) 
8 

the configuration %=0 becomes unstable and decays to its ground state \% \ = a, yielding a graceful 
exit from inflation. Since in this example the ground state of % is not unique, there is the possibility 
of the formation of topological defects at the end of inflation (see [|], ||, ||] for reviews of topological 



defects in cosmology, and the lectures by Polchinski Q31Q for a discussion of how this scenario arises 
in brane inflation models). 

After the slow-roll conditions break down, the period of inflation ends, and the inflaton begins 
to oscillate around its ground state. Due to couplings of (p to other matter fields, the energy of the 
universe, which at the end of the period of inflation is stored completely in <p, gets transferred to 
the matter fields of the particle physics Standard Model. Initially, the energy transfer was described 



8 Note that the slow-roll conditions ( JZ. 1 Q ) and ( 2.1 1 ) were derived assuming that H is given by the contribution of 
<p to V which is not the case here. 
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perturbatively [32, 33]. Later, it was realized [34, 35, 36, 37] that through a parametric resonance 
instability, particles are very rapidly produced, leading to a fast energy transfer ("preheating"). The 
quanta later thermalize, and thereafter the universe evolves as described by SBB cosmology. 



3. Theory of Cosmological Perturbations: A Short Review 

After this review of inflationary cosmology (see e.g. [ |3"8| ] for a more complete recent review), 
we turn to the discussion of the main success of inflationary cosmology, namely the fact that it pro- 
vides a causal mechanism for generating small amplitude inhomogeneities. The reader is referred 



to [ p9[ ] for a comprehensive analysis of this theory of cosmological perturbations, and [ }40[ ] for a 
pedagogical overview, from which most of this section is drawn. 

First, we describe the Newtonian theory of cosmological perturbations, mainly to develop in- 
tuition for the main physical effects. The range of validity of the Newtonian analysis is restricted 
to sub-Hubble scales at late times. In the second subsection, we then summarize the general rel- 
ativistic theory of fluctuations. Linear fluctuations can be trivially quantized. The quantization of 
cosmological perturbations is presented in Subsection 3. 

3.1 Newtonian Theory of Cosmological Perturbations 

The growth of density fluctuations is a consequence of the purely attractive nature of the grav- 
itational force. Given an density excess <5p localized about some point x in space. This inhomo- 
geneity produces a force which attracts the surrounding matter towards x. The magnitude of this 
force is proportional to 8p. Hence, in a non-expanding background, by Newton's second law 

Sp~G5p, (3.1) 

where G is Newton's gravitational constant. Thus, an exponential growth of the fluctuations is 
induced. 

If, as required by consistency in General Relativity, we consider density fluctuations in an 



expanding background, then the expansion of space leads to a friction term in (|3.1[). Hence, instead 
of an exponential instability to the development of fluctuations, the growth of fluctuations will be 
as a power of time. The main goal of the theory of cosmological perturbations is to determine 
how the power-law instability depends on the background cosmology and on the length scale of the 
fluctuations. 



3.2 Perturbations about Minkowski Space-Time 

To develop some physical intuition, we first consider the evolution of hydrodynamical matter 
fluctuations in a fixed non-expanding background. 

In this context, matter is described by a perfect fluid, and gravity by the Newtonian gravita- 
tional potential (p. The fluid variables are the energy density p, the pressure p, the fluid velocity v, 
and the entropy density 5. The basic hydrodynamical equations are 

p + V(pv) = 
v+(v-V p )v + -V p p + V p p = 
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V 2 p (p = AnGp (3.2) 
S+(v-V p )S = 

p = p(p,s), 

where the subscript p indicates that physical as opposed to comoving coordinates are used. The first 
equation is the continuity equation, the second is the Euler (force) equation, the third is the Poisson 
equation of Newtonian gravity, the fourth expresses entropy conservation, and the last describes 
the equation of state of matter. The derivative with respect to time is denoted by an over-dot. 

The background is given by the background energy density p () , the background pressure po, 
vanishing velocity, constant gravitational potential (po and constant entropy density So. Note that 
the background Poisson equation is not satisfied. 

The equations for cosmological perturbations are obtained by perturbing the fluid variables 
about the background, 

p = po + 8p 
v = <5v 

p = po + 8p (3.3) 
<p = <po + S<p 
S = S + 8S, 

where the fluctuating fields Sp ,8\,8p,8(p and SS are functions of space and time, by inserting 



these expressions into the basic hydrodynamical equations (3.2), and by linearizing. After com- 
bining the resulting first order equations, we get the following differential equations for the energy 
density fluctuation <5p and the entropy perturbation 8S 

Sp - c 2 s V 2 p 8p - 4nGp 8p = aV 2 p 8S (3.4) 
SS = 0, 

where the variables c 2 and a describe the equation of state 

8p = c 2 8p + o8S (3.5) 

with 

<* = %X <3 -« 

denoting the square of the speed of sound. 

Since the equations are linear, we can work in Fourier space. Each Fourier component Spk(t) 
of the fluctuation field 8p(x,t) 

Sp(x,t) = [e**8p k (t) (3.7) 



evolves independently. 

The fluctuations can be classified as follows: If 8S vanishes, we have adiabatic fluctuations. 
If the SS is non-vanishing but 8p = 0, we speak of an entropy fluctuation. 



The first conclusions we can draw from the basic perturbation equations (3.4) are that 
1) entropy fluctuations do not grow, 
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2) adiabatic fluctuations are time-dependent, and 

3) entropy fluctuations seed an adiabatic mode. 

Taking a closer look at the equation of motion for 8p, we see that the third term on the left 
hand side represents the force due to gravity, a purely attractive force yielding an instability of flat 



space-time to the development of density fluctuations (as discussed earlier, see (3.1)). The second 
term on the left hand side of (^4j) represents a force due to the fluid pressure which tends to set 
up pressure waves. In the absence of entropy fluctuations, the evolution of Sp is governed by the 
combined action of both pressure and gravitational forces. 



Restricting our attention to adiabatic fluctuations, we see from (3.4) that there is a critical 



wavelength, the Jeans length, whose wavenumber kj is (in physical coordinates) given by 

kj=(^) 1/2 . (3-8) 
Fluctuations with wavelength longer than the Jeans length (k <^kj) grow exponentially 

Sp k (t) ~ e mk ' with co k ~ 4(7iGp ) 1/2 (3.9) 

whereas short wavelength modes (k S> kj) oscillate with frequency (Ou ~ c s k. Note that the value of 
the Jeans length depends on the equation of state of the background. For a background dominated 
by relativistic radiation, the Jeans length is large (of the order of the Hubble radius H~ l {t)), whereas 
for pressure-less matter it goes to zero. 

Next, we study Newtonian cosmological fluctuations about an expanding background. In this 
case, the background equations are consistent (the non-vanishing average energy density leads to 
cosmological expansion). However, we are neglecting general relativistic effects (the fluctuations 
of the metric) which dominante on length scales larger than the Hubble radius H~ l (t). 

The background cosmological model is given by the energy density Po(f), the pressure po(t), 
and the recessional velocity vo = H(t)x p where x p is the physical coordinate vector. The space- 
and time-dependent fluctuating fields are defined in analogy to the previous section: 

p(t,x) = po(t)(l + 8 e (t,x)) 

v(t,x) = y (t,x) + 8\{t,x) (3.10) 
p(t,x) = p (t) + 8p{t,x), 

where 8 e is the fractional energy density perturbation (we are interested in the fractional rather than 
in the absolute energy density fluctuation!), and the pressure perturbation 8p is defined as in ([331) . 
In addition, there is the possibility of a non- vanishing entropy perturbation defined as in (EO). 



We now insert this ansatz into the basic hydrodynamical equations ( |3.2| ), linearize in the per- 
turbation variables, and combine the first order differential equations for 8 e and 8p into a single 
second order differential equation for 8 e . The result simplifies if we work in comoving coordinates 
x. After some algebra, we obtain the following equation which describes the time evolution of 
density fluctuations: 

2 

8 e + 2H8 e -^V 2 8 e -47iGp 8 e = -%8S, (3.11) 
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where V is the partial derivative vector with respect to comoving coordinates. In addition, we have 
the equation of entropy conservation 

55 = 0. (3.12) 



Comparing with the equations (3.4) obtained in the absence of an expanding background, 
we see that the only difference is the presence of a Hubble damping term in the equation for 8 e . 
This term will moderate the exponential instability of the background to long wavelength density 
fluctuations. In addition, it will lead to a damping of the oscillating solutions on short wavelengths. 



More specifically, for physical wavenumbers k p <C kj (where kj is again given by (|3.8|)), and in 



a matter-dominated background cosmology, the general solution of ( p.ll| ) in the absence of any 
entropy fluctuations is given by 

4(0 =cif 2/3 + c 2 r 1 , (3.13) 

where and c 2 are two constants determined by the initial conditions, and we have dropped the 
subscript £ in expressions involving 8 e . There are two fundamental solutions, the first a growing 
mode with 4(0 ~ a(t), the second a decaying mode with 4(0 ~ f -1 . On short wavelength, one 
obtains damped oscillatory motion: 

4(0 ~ a- 1/2 (t)exp(±ic s k f dt'ar x (t')) . (3.14) 

Before going on to the relativitic theory of cosmological perturbations, we will pause to intro- 
duce terminology used in cosmology to describe the fluctuations. 

Let us consider perturbations on a fixed comoving length scale given by a comoving wavenum- 
ber k. The corresponding physical length increases as a{t). This is to be compared to the Hubble 
radius H~ l (t) which scales as t provided a{t) grows as a power of t. In the late time Universe, 
a(t) ~ f 1 / 2 in the radiation-dominated phase (i.e. for t < t eq ), and a(t) ~ t 2 ^ in the matter- 
dominated period (t eq < t < t ). Thus, at sufficiently early times, all comoving scales had a physical 
length larger than the Hubble radius. If we consider large cosmological scales (e.g. those corre- 
sponding to the observed CMB anisotropics or to galaxy clusters), the time tn(k) of "Hubble radius 
crossing" (when the physical length was equal to the Hubble radius) was in fact later than t eq . The 
time of Hubble radius crossing plays an important role in the evolution of cosmological perturba- 
tions. 

Cosmological fluctuations can be described either in position space or in momentum space. 
In position space, we compute the root mean square mass fluctuation 8M /M(k,t) in a sphere of 
radius I = 2%jk at time t. A scale-invariant spectrum of fluctuations is defined by the relation 

^-(k,t H (k)) = const. . (3.15) 

Such a spectrum was first suggested by Harrison [ pd| ] and Zeldovich [12] as a reasonable choice 
for the spectrum of cosmological fluctuations. The "spectral index" n of cosmological fluctuations 
is defiined by the relation 

(^-) 2 (k,t H (k))-k"-K (3.16) 
Thus, a scale-invariant spectrum corresponds to n = 1 . 
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To go to momentum space representation, the fractional spatial density contrast is expanded in 
a Fourier series: 



Se(x,t) 



d 3 k§AKt)e ikx . 



The power spectrum P§ of density fluctuations is defined by 

P s (k) = ^\8 e (k)\ 2 , 



(3-17) 



(3.18) 



where k is the magnitude of k. For simplicity, the distribution of fluctuations is taken to be Gaussian 
so that the fluctuation amplitude only depends on k. 

The power spectrum of the gravitational potential (p is defined by 



P< P (k) = ^\S~(p(k)\ 2 . 
The two power spectra ( 3.18 ) and ( 3.19| ) are related by the Poisson equation (3.2) 

P v (k) ~ k- 4 P 5 (k). 



(3.19) 



(3.20) 



The condition of scale-invariance can be expressed in terms of the power spectrum evaluated 
at a fixed time. To obtain this condition, we first use the time dependence of the fractional density 
fluctuation from ( 3.13 ) to determine the mass fluctuations at a fixed time. We need to use the fact 
that the time of Hubble radius crossing is given by 



a(t H (k))k = pt H (k), 



(3.21) 



where /3 = 2 or /3 = 3/2 in the radiation and matter dominated phases, respectively. Making use of 
(06b we find 

(3.22) 



n+3 



(—Y(k,t) ~ k 

Since, for reasonable values of the index of the power spectrum, SM/M(k,t) is dominated by the 
Fourier modes with wavenumber k, we find that ( 3.22 ) implies 

P S (k)~k" +3 . 

or, equivalently, 



(3.23) 



rJl-1 



(3.24) 



P v {k) ~ *" 

3.3 Relativistic Theory of Cosmological Fluctuations 

The Newtonian theory of cosmological fluctuations discussed in the previous section breaks 
down on scales larger than the Hubble radius because it neglects perturbations of the metric, and 
because on large scales the metric fluctuations dominate the dynamics. 

To show why metric fluctuations are important on scales larger than the Hubble radius, we 
can use a "separate universe" argument. On such scales, one should be able to approximately de- 
scribe the evolution of the space-time by applying the first FRW equation (2.1) of homogeneous 
and isotropic cosmology to the local Universe (this approximation is made more rigorous in [43]). 
Based on this equation, a large-scale fluctuation of the energy density will lead to a fluctuation 
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("8a") of the scale factor a which grows in time, a manifestation of the self gravitational amplifi- 
cation of fluctuations on length scales A greater than the Hubble radius. 

Let us now turn to the rigorous analysis of cosmological fluctuations in the context of general 
relativity, where both metric and matter inhomogeneities are taken into account. We will consider 



fluctuations about a homogeneous and isotropic background cosmology, given by the metric (|1.3|), 
which can be written in conformal time tj (defined by dt = a(t)dr\) as 

ds 2 =a{ri) 2 (dri 2 -dx 2 ). (3.25) 

The theory of cosmological perturbations is based on expanding the Einstein equations to lin- 
ear order about the background metric. The theory was initially developed in pioneering works 



by Lifshitz [44]. Significant progress in the understanding of the physics of cosmological fluctu- 



ations was achieved by Bardeen [45] who realized the importance of subtracting gauge artifacts 



(see below) from the analysis (see also [46]). The following discussion is based on Part I of the 
comprehensive review article [39]. Other reviews - emphasizing different aspects or approaches - 
are @f|| ff9[ [50|]. 

The first step in the analysis of metric fluctuations is to classify them according to their trans- 
formation properties under spatial rotations. There are scalar, vector and second rank tensor fluc- 
tuations. In linear theory, there is no coupling between the different fluctuation modes, and hence 
they evolve independently (for some subtleties in this classification, see [|5l|]). 



We begin by expanding the metric about the FRW background metric g^l given by ( |3.25| ): 



8nv = sfl + s 8»v ■ (3.26) 

The background metric depends only on time, whereas the metric fluctuations 8g^ v depend on 
both space and time. Since the metric is a symmetric tensor, there are at first sight 10 fluctuating 
degrees of freedom in 8g^ v . 

There are four degrees of freedom which correspond to scalar metric fluctuations (the only 
four ways of constructing a metric from scalar functions): 

where the four fluctuating degrees of freedom are denoted (following the notation of [[3^]) (j),B,E, 
and \\f, a comma denotes the ordinary partial derivative (if we had included spatial curvature of 
the background metric, it would have been the covariant derivative with respect to the background 
spatial metric), and 8{j is the Kronecker symbol. 

There are four vector degrees of freedom of metric fluctuations, consisting of the four ways of 
constructing metric fluctuations from three vectors: 




where S; and Fj are two divergence-less vectors (for a vector with non-vanishing divergence, the 
divergence contributes to the scalar gravitational fluctuation modes). 



15 



Topics in Cosmology 



Robert Brandenberger 



Finally, there are two tensor modes which correspond to the two polarization states of gravita- 
tional waves: 

Sg^ v = -a 2 ( ° n ,° | , (3.29) 



Ohij 



where h{j is trace-free and divergence-less 



h\ = hjj = 0. (3.30) 

Gravitational waves do not couple at linear order to the matter fluctuations. Vector fluctuations 
decay in an expanding background cosmology and hence are not usually cosmologically important. 
Thus, the most important fluctuations, at least in inflationary cosmology, are the scalar metric 
fluctuations, the fluctuations which couple to matter inhomogeneities and which are the relativistic 
generalization of the Newtonian perturbations considered in the previous section. 

The theory of cosmological perturbations is at first sight complicated by the issue of gauge in- 
variance. The coordinates t, x of space-time carry no independent physical meaning. By performing 
a small-amplitude transformation of the space-time coordinates (called "gauge transformation" in 
the following), we can create "fictitious" fluctuations in a homogeneous and isotropic Universe. 
These modes are "gauge artifacts". 

In the following we take an "active" view of gauge transformation. Consider two space-time 
manifolds, one of them a homogeneous and isotropic Universe ^#o, the other a physical Universe 
^# with inhomogeneities. A choice of coordinates can be considered as a mapping between the 
manifolds and jtft. A second mapping S will map the same point in j^q into a different point 
in Jt. Using the inverse of these maps S and §}, we can assign two different sets of coordinates 
to points in 

Consider now a physical quantity Q (e.g. the Ricci scalar) on , and the corresponding 
physical quantity Q(°) on j^q Then, in the first coordinate system given by the mapping S, the 
perturbation 8Q of Q at the point p £ ^ is defined by 

SQ(p) = Q(p)-Q {0) {S-\p)). (3.31) 

In the second coordinate system given by S, the perturbation is defined by 

8~Q(p) =Q{p)-Q {0) {®-\p)). (3.32) 

The difference 

AQ(p) = 8~Q(p) - 8Q(p) (3.33) 

is a gauge artifact and carries no physical significance. 

Some of the metric perturbation degrees of freedom introduced in the first subsection are 
thus gauge artifacts. To isolate these, we must study how coordinate transformations act on the 
metric. There are four independent gauge degrees of freedom corresponding to the coordinate 
transformations 

+ (3.34) 

The time component of ^ leads to a scalar metric fluctuation. The spatial three vector E, 1 can 
be decomposed as 

= & + (3-35) 
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(where Y J is the spatial background metric) into a transverse piece Q r which has two degrees of 
freedom which yield vector perturbations, and the second term (given by the gradient of a scalar t, ) 
which gives a scalar fluctuation. Thus, there are two scalar gauge modes given by ^° and and 
two vector modes given by the transverse three vector £ ( ' r Thus, there remain two physical scalar 
and two vector fluctuation modes. The gravitational waves are gauge-invariant. 

Let us now focus on how the scalar gauge transformations (i.e. the transformations given by 
(§° and act on the scalar metric fluctuation variables (j),B,E, and y. An immediate calculation 
yields: 

; = 0-^°-(£ )' 

a 

B = B + %°-£ (3.36) 
E = E-% 

a ' f o 

¥ = V+-<T, 
a 

where a prime indicates the derivative with respect to conformal time tj . 

There are two approaches to deal with the gauge ambiguities. The first is to fix a gauge, i.e. to 
pick conditions on the coordinates which completely eliminate the gauge freedom, the second is to 
work with a basis of gauge-invariant variables. 

If one wants to adopt the gauge-fixed approach, there are many different gauge choices. Note 
that the often used synchronous gauge determined by 8g 0fl = does not totally fix the gauge. A 
convenient system which completely fixes the coordinates is the so-called longitudinal or confor- 
mal Newtonian gauge defined by B = E = 0. 

If one prefers a gauge-invariant approach, there are many choices of gauge-invariant variables. 
A convenient basis first introduced by [ p3| ] is the basis given by 

O = (j) + -[(B-E')a]' (3.37) 

y = W --(B-E'). (3.38) 
a 

The gauge-invariant variables <I> and *P coincide with the corresponding diagonal metric perturba- 
tions <p and y m longitudinal gauge. 



Note that the variables defined in (3.37) are gauge-invariant only under linear space-time co- 



ordinate transformations. Beyond linear order, the structure of perturbation theory becomes much 



more involved. In fact, one can show [ ]52| ] that the only fluctuation variables which are invariant 
under all coordinate transformations are perturbations of variables which are constant in the back- 
ground space-time. Beyond linear order there is also mixing of scalar, vector and tensor modes. 

To derive the equations of motion for the fluctuations, the starting point is the set of Einstein 
equations 

G^v = SkGT^ v , (3.39) 

where G^ v is the Einstein tensor associated with the space-time metric g^ v , and 7^ v is the energy- 
momentum tensor of matter. We insert the ansatz for metric and matter perturbed about a FRW 
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background 9o(v))- 

gfivfan) = 8 { !1(V) + Sg^v(x,v) (3-40) 
<p(x,Tj) = <po(rj) + S<p(x,rj), (3.41) 

(where we have for simplicity replaced general matter by a scalar matter field q>) and expand to 
linear order in the fluctuating fields, obtaining the following equations: 

<5G MV = SnGdT^ . (3.42) 

In the above, 8g^ v is the perturbation in the metric and 8(p is the fluctuation of the matter field (p. 
In a gauge-fixed approach, one can start with the metric in longitudinal gauge 

ds 2 = a 2 [(l+2<j))dri 2 - (l-2\if)y i jdjc i dx j ] (3.43) 

and insert this ansatz into the general perturbation equations ( |3.42| ). This yields the following set 
of equations of motion: 

-3Jf(jr0 + y/) +V 2 yf = 4nGa 2 8T$ 

(Jt?<p + y) i = 4nGa 2 8T, Q (3.44) 

[ (2 jr' + jr 2 ) + 4 + y" + 2 y] sj 

+^ Dd i_l y ik Dk , = -tocJsTj, 

where D = (j) — y/ and Jif = a' /a. 

If no anisotropic stress is present in the matter at linear order in fluctuating fields, i.e. if 
8Tj = for i ^ j, then the two metric fluctuation variables coincide, i.e. <p = y . This will be the 
case in most simple cosmological models, e.g. in theories with matter described by a set of scalar 
fields with canonical form of the action, and in the case of a perfect fluid with no anisotropic stress. 

In the simple case in which matter described in terms of a single scalar field (p, then in longi- 
tudinal gauge ( 3.44 ) reduce to the following set of equations of motion 

V 2 ^)-3J^<p' - (jf' + 2^ 2 )0 = 4xG(%8(p' + Va 2 8(p) 

4 + = 4nG(p 8(p (3.45) 
tp" + 3J^f<p' + + 2Jf 2 )<p = 47iG(%8(p' -Va 2 8(p), 

where V denotes the derivative of V with respect to (p. These equations can be combined to give 
the following second order differential equation for the relativistic potential : 

</> +2 I ..tf - — \(? -V 2 + 2 I Jf'-jr% I = 0. (3.46) 




This is the final result for the classical evolution of cosmological fluctuations 9 . 



9 Note that we have implicitly assumed that the background matter field is slowly rolling, as it does in slow-roll 
inflation. In the case that it is time-independent, then the leading metric fluctuations are quadratic in the matter inhomo- 
geneities, as discussed in [1531] . 
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There are similarities between the above equation of motion for the relativistic perturbations 



and the equation ( |3.11| ) obtained in the Newtonian theory. The final term in ( p.46| ) is the force 
due to gravity leading to the instability, the second to last term is the pressure force leading to 
oscillations (relativistic since we are considering matter to be a relativistic field), and the second 
term is the Hubble friction term. For each wavenumber there are two fundamental solutions. On 
small scales (k > H), the solutions correspond to damped oscillations, on large scales (k < H) the 
oscillations freeze out and the dynamics is governed by the gravitational force competing with the 
Hubble friction term. Note, in particular, how the Hubble radius naturally emerges as the scale 
where the nature of the fluctuating modes changes from oscillatory to frozen. 

Considering the equation in a bit more detail, observe that if the equation of state of the back- 
ground is independent of time (which will be the case if M 3 =%= 0), then in an expanding 



background, the dominant mode of ( |3.4q ) is constant, and the sub-dominant mode decays. If the 
equation of state is not constant, then the dominant mode is not constant in time. Specifically, at 
the end of inflation Jrif < 0, and this leads to a growth of <j) (see below). 



To study the quantitative implications of the equation of motion (3.46), it is convenient to in- 
troduce [54, 55] the variable £ (which, up to correction term of the order V 2 which is unimportant 
for large-scale fluctuations, is equal to the curvature perturbation M in comoving gauge [p§]) by 



where 



3 l+w 



P 

w = — 
9 



(3.47) 



(3.48) 



characterizes the equation of state of matter. In terms of £, the equation of motion (3.46) takes on 
the form 



3 4 



(3.49) 



On large scales, the right hand side of the equation is negligible, which leads to the conclusion that 
large-scale cosmological fluctuations satisfy 



£(l+w) = 0. 



(3.50) 



This implies that £ is constant except possibly if 1 + w = at some point in time during the 
cosmological evolution (which occurs during reheating in inflationary cosmology if the inflaton 
field undergoes oscillations - see [|7|] and [58], 59] for discussions of the consequences in single and 
double field inflationary models, respectively). In single matter field models it is indeed possible 
to show that £ = on super-Hubble scales independent of assumptions on the equation of state 



[ ]60| , |61Q . This "conservation law" makes it easy to relate initial fluctuations to final fluctuations in 
inflationary cosmology, as will be illustrated in the following. 

Consider an application to inflationary cosmology. We return to the space-time sketch of 
the evolution of fluctuations - see Figure ([]]) - and use the conservation law ( 3.50D - in the form 
£ = const on large scales - to relate the amplitude of at initial Hubble radius crossing during the 
inflationary phase (at t = t{(k)) with the amplitude at final Hubble radius crossing at late times (at 
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t = tf(k)). Since both at early times and at late times = on super-Hubble scales, as the equation 
of state is not changing, ( 3.50| ) and ( 3.47 ) lead to 



If the initial values of the perturbations are known, then the above equation allows us to evaluate 
the fluctuation amplitude at the time tf(k) of re-entry into the Hubble radius. 



The time-time perturbed Einstein equation (the first equation of ( 3.44 )) relates the value of (j) 
at the initial Hubble radius crossing to the amplitude of the fractional energy density fluctuations at 
that time. This, together with the fact that the amplitude of the scalar matter field quantum vacuum 
fluctuations is of the order H, yields 

V 

${ti{k))~H—(tt{k)). (3.52) 
In the late time radiation dominated phase, w = 1/3, whereas during slow-roll inflation 

l+w(f,-(*))=;^(ff(*))- (3.53) 

Making, in addition, use of the slow roll conditions satisfied during the inflationary period 

3H<po ~ -V 
, 87TG 

H 2 ~ — V, (3.54) 

v 3 / 2 

<t>(t f (k)) ~ =p-3-fe(*)), (3.55) 
Vm P i 

which gives the position space amplitude of cosmological fluctuations on a scale labelled by the 
comoving wavenumber k at the time when the scale re-enters the Hubble radius at late times, a 



we arrive at the final result 



result first obtained in the case of the Starobinsky model [15] of inflation in [16], and later in the 



context of scalar field-driven inflation in 062], p3L p4L 54]. 



In the case of slow roll inflation, the right hand side of ( 3.55p is, to a first approximation, 



independent of k, and hence the resulting spectrum of fluctuations is nearly scale-invariant. 

3.4 Quantum Theory of Cosmological Fluctuations 

In many models of the very early Universe, in particular in inflationary cosmology, primordial 
inhomogeneities emerge from quantum vacuum fluctuations on microscopic scales (wavelengths 
smaller than the Hubble radius). The wavelength is then stretched relative to the Hubble radius, 
becomes larger than the Hubble radius at some time and the perturbation then propagates on super- 
Hubble scales until re-entering at late cosmological times. In the context of a Universe with a de 



Sitter phase, the quantum origin of cosmological fluctuations was first discussed in Jig]. In an in- 
flationary universe, it is easy to justify focusing attention on the quantum fluctuations: any classical 
fluctuations present at the beginning of inflation are red-shifted during the period of inflation, and 
will thus be irrelevant to scales probed in observations today. On small scales, a quantum vacuum 
remains. 
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To understand the role of the Hubble radius, consider the equation of a free scalar matter field 
<p on an unperturbed expanding background: 

V 2 

<i> + m<p-^r<p = 0. (3.56) 

a 1 

The second term on the left hand side of this equation leads to damping of cp with a characteristic 
decay rate given by H. As a consequence, in the absence of the spatial gradient term, <p would be 
of the order of magnitude Hq>. Thus, comparing the second and the third terms on the left hand 
side, we immediately see that the microscopic (spatial gradient) term dominates on length scales 
smaller than the Hubble radius, leading to oscillatory motion, whereas this term is negligible on 
scales larger than the Hubble radius, and the evolution of q> is determined primarily by gravity. 
Note that in general cosmological models the Hubble radius is much smaller than the horizon (the 
forward light cone calculated from the initial time). In an inflationary universe, the horizon is 
larger by a factor of exp(Af), where N is the number of e-foldings of inflation. It is very important 
to realize this difference, a difference which is obscured in most articles on cosmology in which 
the term "horizon" is used when "Hubble radius" is meant. Note, in particular, that the homoge- 
neous inflaton field contains causal information on super-Hubble but sub-horizon scales. Hence, 



it is completely consistent with causality Q57| ] to have a microphysical process related to the back- 
ground scalar matter field lead to exponential amplification of the amplitude of fluctuations during 
reheating on such scales, as it does in models in which entropy perturbations are present and not 



suppressed during inflation [p8|, |59|]. 

To understand the generation and evolution of fluctuations in current models of the very early 
Universe, we need both Quantum Mechanics and General Relativity, i.e. quantum gravity. At first 
sight, we are thus faced with an intractable problem, since the theory of quantum gravity is not 
yet established. We are saved by the fact that today on large cosmological scales the fractional 
amplitude of the fluctuations is smaller than 1. Since gravity is a purely attractive force, the fluctu- 
ations had to have been - at least in the context of an eternally expanding background cosmology 
- very small in the early Universe. Thus, a linearized analysis of the fluctuations (about a classical 
cosmological background) is self-consistent. 

From the classical theory of cosmological perturbations discussed in the previous subsection, 
it follows that the analysis of scalar metric inhomogeneities can be reduced - after extracting gauge 
artifacts - to the study of the evolution of a single fluctuating variable. Thus, the quantum the- 
ory of cosmological perturbations must be reducible to the quantum theory of a single free scalar 
field which we will denote by v. Since the background in which this scalar field evolves is time- 
dependent, the mass of v will be time-dependent. The time-dependence of the mass will lead to 
quantum particle production over time if we start the evolution in the vacuum state for v. As we will 
see, this quantum particle production corresponds to the development and growth of the cosmologi- 
cal fluctuations. The quantum theory of cosmological fluctuations provides a consistent framework 
to study both the generation and the evolution of metric perturbations. 

In order to obtain the action for linearized cosmological perturbations, we expand the action 
to quadratic order in the fluctuating degrees of freedom. The linear terms cancel because the 
background is taken to satisfy the background equations of motion. 
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We begin with the Einstein-Hilbert action for gravity and the action of a scalar matter field 

S = J d A x^[-j^R+ l -d^(p-V((p)], (3.57) 

where R is the Ricci curvature scalar. 

The simplest way to proceed is to work in longitudinal gauge. The next step is to reduce 
the number of degrees of freedom. Since for scalar field matter there are no anisotropic stresses to 
linear order, the off-diagonal spatial Einstein equations force y=<$>. The two remaining fluctuating 
variables and (p are linked by the Einstein constraint equations since there cannot be matter 
fluctuations without induced metric fluctuations. 



The two nontrivial tasks of the lengthy \p% computation of the quadratic piece of the action 
is to find out what combination of <p and yields the variable v in terms of which the action has 
canonical kinetic term, and what the form of the time-dependent mass is. In the context of scalar 
field matter, the quantum theory of cosmological fluctuations was developed by Mukhanov [65, 66] 
and Sasaki [|67|]. The result is the following form of the action quadratic in the perturbations: 

S (2) = lfd 4 x[v' 2 -v. i v, i + Z jv 2 ], (3.58) 



a[8(p + ^(j)] (3.59) 

(3.60) 



where the canonical variable v (the "Sasaki-Mukhanov variable" introduced in [ ]66| ] - see also [p8|]) 
is given by 

% 

Jif' 

and 

= ^Po 

Z J? 

As long as the equation of state does not change over time and <p are proportional and 
hence 

z{r\)~a{r\). (3.61) 

Note that the variable v is related to the curvature perturbation M in comoving coordinates intro- 
duced in [56] and closely related to the variable £ used in |55[l: 

v = Z M. (3.62) 



The equation of motion which follows from the action ( p.58[ ) is (in momentum space) 

// 

v" k + k 2 v k --v k = 0, (3.63) 

z 



where v k is the k'th Fourier mode of v. As a consequence of ( p.61| ), the tachyonic mass term in the 
above equation is given by the Hubble scale 

// 

k 2 H = -~H 2 . (3.64) 



Thus, it immediately follows from (3.63) that on small length scales, i.e. for k > ku, the solutions 
for Vk are constant amplitude oscillations . These oscillations freeze out at Hubble radius crossing, 
i.e. when k = k^. On longer scales (k ^ %), the solutions for v k increase as z: 

Vk ~ z , k <C kn ■ (3.65) 
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The state of the fluctuations becomes a squeezed quantum state. 



Given the action ( 3.58 ), the quantization of the cosmological perturbations can be performed 
by canonical quantization (in the same way that a scalar matter field on a fixed cosmological back- 
ground is quantized [|6^]). 

The final step in the quantum theory of cosmological perturbations is to specify an initial 
state. Since in inflationary cosmology all pre-existing classical fluctuations are red-shifted by the 
accelerated expansion of space, one usually assumes (we will return to a criticism of this point 
when discussing the trans-Planckian problem of inflationary cosmology) that the field v starts out 
at the initial time t{ mode by mode in its vacuum state. Two questions immediately emerge: what 
is the initial time t(, and which of the many possible vacuum states should be chosen. It is usually 
assumed that since the fluctuations only oscillate on sub-Hubble scales, the choice of the initial 
time is not important, as long as it is earlier than the time when scales of cosmological interest 
today cross the Hubble radius during the inflationary phase. The state is usually taken to be the 
Bunch-Davies vacuum (see e.g. [|6^]), since this state is empty of particles at t,- in the coordinate 



frame determined by the FRW coordinates (see e.g. [J70p for a discussion of this point), and since 
the Bunch-Davies state is a local attractor in the space of initial states in an expanding background 
(see e.g. [f7l|]). Thus, we choose the initial conditions 

v k (r]i) = (3.66) 



vM) = 

where here (Ok = k, and T] ( is the conformal time corresponding to the physical time 

Let us briefly summarize the quantum theory of cosmological perturbations. In the linearized 
theory, fluctuations are set up at some initial time tj mode by mode in their vacuum state. While the 
wavelength is smaller than the Hubble radius, the state undergoes quantum vacuum fluctuations. 
The accelerated expansion of the background redshifts the length scale beyond the Hubble radius. 
The fluctuations freeze out when the length scale is equal to the Hubble radius. On larger scales, 
the amplitude of increases as the scale factor. This corresponds to the squeezing of the quantum 
state present at Hubble radius crossing (in terms of classical general relativity, it is self-gravity 
which leads to this growth of fluctuations). As discussed e.g. in [f72|], the squeezing of the quantum 
vacuum state sets up the classical correlations in the wave function of the fluctuations which are an 
essential ingredient in the classicalization of the perturbations. Squeezing also leads to the phase 
coherence of fluctuations on all scales, which in turn is responsible for producing the oscillations 
in the CMB angular power spectrum, as realized long before the advent of inflationary cosmology 
©0]- 



m 



We end this subsection with a calculation of the spectrum of curvature fluctuations in infla- 
tionary cosmology. 

We need to compute the power spectrum &,$>(k) of the curvature fluctuation & defined in 



(3.62). The idea in calculating the power spectrum at a late time t is to first relate the power 



spectrum via the growth rate (3.65) of v on super-Hubble scales to the power spectrum at the time 



tn(k) of Hubble radius crossing, and to then use the constancy of the amplitude of v on sub-Hubble 
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scales to relate it to the initial conditions (3 

.2 



). Thus 



2r 



,0 = ew(t) = *v z (*)|v*(0 



(3.67) 



= ^- a w(^y)W(*))P 

= k\- 2 (t H (k))\n(t H (k))\ 2 

~ fcV 2 (fe(£))|v*(* ; )| 2 , 

where in the final step we have used ( |3.61| ) and the constancy of the amplitude of v on sub-Hubble 
scales. Making use of the condition a _1 (f#(A:))fc = H for Hubble radius crossing, and of the initial 
conditions ( 3.66 ), we immediately see that 



(3.68) 



and that thus a scale invariant power spectrum with amplitude proportional to H 2 results. 

The quantization of gravitational waves parallels the quantization of scalar metric fluctuations, 
but is more simple because there are no gauge ambiguities. The starting point is the action ( 3.57 ), 
into which we insert the metric which corresponds to a classical cosmological background plus 
tensor metric fluctuations: 



ds 2 = a 2 (r\) \dt] 2 — (8ij + h i j)dx t dx-'] , 



(3.69) 



where the second rank tensor /z, ; (t],x) represents the gravitational waves, and in turn can be de- 
composed as 



hij(r),x) = /z + (7],x)4 + ^(7],x)ef. 



(3.70) 



'IJ I ■~A.\-ll-/-lJ 

into the two polarization states. Here, and are two fixed polarization tensors, and h + and h x 
are the two coefficient functions. 

To quadratic order in the fluctuating fields, the action consists of separate terms involving h + 
and h r . Each term is of the form 



S (2) = Jd^[h> 2 -(Vh) 2 }, 



leading to the equation of motion 



a' . / 



h k + 2-h k + k z h k = 0. 



The variable in terms of which the action (3.71) has canonical kinetic term is 



and its equation of motion is 



JU* = ah . 



H k + {k 2 -—)n k = 0. 



(3.71) 



(3.72) 



(3.73) 



(3.74) 



This equation is very similar to the corresponding equation (3.63) for scalar gravitational inhomo- 
geneities, except that in the mass term the scale factor a(rj) replaces z(tj), which leads to a very 
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different evolution of scalar and tensor modes during the reheating phase in inflationary cosmology 
during which the equation of state of the background matter changes dramatically 

Based on the above discussion we have the following theory for the generation and evolution 
of gravitational waves in an accelerating Universe (first developed by Grishchuk [|75|]): waves exist 
as quantum vacuum fluctuations at the initial time on all scales. They oscillate until the length 
scale crosses the Hubble radius. At that point, the oscillations freeze out and the quantum state of 
gravitational waves begins to be squeezed in the sense that 

/*t(i7) (3-75) 



which, from ( |3.73| ) corresponds to constant amplitude of The squeezing of the vacuum state 
leads to the emergence of classical properties of this state, as in the case of scalar metric fluctua- 
tions. 



4. Towards a New Paradigm: String Gas Cosmology 

In spite of the phenomenological successes of inflationary cosmology in addressing some key 
problems of standard cosmology and in providing a predictive theory of structure formation, current 
models of inflation are faced with some key conceptual problems which motivate the search for a 
new early universe paradigm based on new fundamental physics, string theory being the most 
promising candidate. Below, we first list some of the key problems for inflation, and then discuss a 
toy model designed to explore possible cosmological consequences of some of the key new features 
of string theory. 



4.1 Conceptual Problems of Scalar Field-Driven Inflation 

Nature of the Inflaton 

In the context of General Relativity as the theory of space-time, matter with an equation of 
state p ~ — p is required in order to obtain almost exponential expansion of space. If we describe 
matter in terms of fields with canonical kinetic terms, a scalar field is required since in the context 
of usual renormalizable field theories it is only for scalar fields that a potential energy function in 
the Lagrangian is allowed, and of all energy terms only the potential energy can yield the required 
equation of state. 

In order for scalar fields to generate a period of cosmological inflation, the potential energy 
needs to dominate over the kinetic and spatial gradient energies. It is generally assumed that spatial 
gradient terms can be neglected. This is, however, not true in general. Next, assuming a homo- 
geneous field configuration, we must ensure that the potential energy dominates over the kinetic 
energy. This leads to the first "slow-roll" condition. Requiring the period of inflation to last suf- 
ficiently long leads to a second slow-roll condition, namely that the (p term in the Klein-Gordon 
equation for the inflaton (p be negligible. Scalar fields charged with respect to the Standard Model 
symmetry groups do not satisfy the slow-roll conditions. 

Assuming that both slow-roll conditions hold, one obtains a "slow-roll trajectory" in the phase 
space of homogeneous cp configurations. In large-field inflation models such as "chaotic inflation" 
[27] and "hybrid inflation" [RO], the slow -roll trajectory is a local attractor in initial condition space 
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[19] (even when linearized metric perturbations are taken into account [29]), whereas this is not 



the case [ ]2Q[ | in small-field models such as "new inflation" [|24|]. As shown in [|26|], this leads to 
problems for some models of inflation which have recently been proposed in the context of string 
theory. To address this problem, it has been proposed that inflation may be future-eternal [76] and 
that it is hence sufficient that there be some configurations in initial condition space which give rise 
to inflation within one Hubble patch, inflation being then self-sustaining into the future. However, 
one must still ensure that slow-roll inflation can locally be satisfied. 

Many models of particle physics beyond the Standard Model contain a plethora of new scalar 
fields. One of the most conservative extensions of the Standard Model is the MSSM, the "Minimal 
Supersymmetric Standard Model". According to a recent study, among the many scalar fields in 
this model, only a hand-full can be candidates for a slow-roll inflaton, and even then very special 
initial conditions are required [f77[]. The situation in supergravity and superstring-inspired field 
theories may be more optimistic, but the issues are not settled (see e.g. [Q, [| ^ for recent reviews). 

Hierarchy Problem 

Assuming for the sake of argument that a successful model of slow-roll inflation has been 
found, one must still build in a hierarchy into the field theory model in order to obtain an accept- 
able amplitude of the density fluctuations (this is sometimes also called the "amplitude problem"). 
Unless this hierarchy is observed, the density fluctuations will be too large and the model is obser- 
vationally ruled out. 

In a wide class of inflationary models, obtaining the correct amplitude requires the introduction 
of a hierarchy in scales, namely [78] 

V(q>) n 

^<10- 12 , (4.1) 

where Acp is the change in the inflaton field during one Hubble expansion time (during inflation), 
and V((p) is the potential energy during inflation. 

This problem should be contrasted with the success of topological defect models (see e.g. 
[Q, |5| [J for reviews) in predicting the right oder of magnitude of density fluctuations without 
introducing a new scale of physics. The GUT scale as the scale of the symmetry breaking phase 
transition (which produces the defects) yields the correct magnitude of the spectrum of density 
fluctuations [79]. Topological defects, however, cannot be the prime mechanism for the origin of 
fluctuations since they do not give rise to coherent adiabatic fluctuations and hence fail to yield 
acoustic oscillations in the angular power spectrum of the CMB anisotropics [80]. 

Trans-Planckian Problem 



A more serious problem is the "trans-Planckian problem" J81Q. Returning to the space-time 
diagram of Figure 1 , we can immediately deduce that, provided that the period of inflation lasted 
sufficiently long (for GUT scale inflation the number is about 70 e-foldings), then all scales inside 
of the Hubble radius today started out with a physical wavelength smaller than the Planck scale at 
the beginning of inflation. Now, the theory of cosmological perturbations is based on Einstein's 
theory of General Relativity coupled to a simple semi-classical description of matter. It is clear 
that these building blocks of the theory are inapplicable on scales comparable and smaller than the 
Planck scale. Thus, the key successful prediction of inflation (the theory of the origin of fluctu- 
ations) is based on suspect calculations since new physics must enter into a correct computation 
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of the spectrum of cosmological perturbations. The key question is as to whether the predictions 
obtained using the current theory are sensitive to the specifics of the unknown theory which takes 
over on small scales. 

One approach to study the sensitivity of the usual predictions of inflationary cosmology to 
the unknown physics on trans-Planckian scales is to study toy models of ultraviolet physics which 
allow explicit calculations. The first approach which was used [82|, 83] is to replace the usual linear 
dispersion relation for the Fourier modes of the fluctuations by a modified dispersion relation, a 
dispersion relation which is linear for physical wavenumbers smaller than the scale of new physics, 
but deviates on larger scales. Such dispersion relations were used previously to test the sensitivity of 
black hole radiation on the unknown physics of the ultra-violet [Q, |^]. It was found [82] that if the 
evolution of modes on the trans-Planckian scales is non-adiabatic, then substantial deviations of the 
spectrum of fluctuations from the usual results are possible. Non-adiabatic evolution turns an initial 
state minimizing the energy density into a state which is excited once the wavelength becomes 
larger than the cutoff scale. Back-reaction effects of these excitations may limit the magnitude of 
the trans-Planckian effects, but - based on our recent study [ 86] - not to the extent initially expected 

From the point of view of fundamental physics, the trans-Planckian problem is not a problem. 
Rather, it yields a window of opportunity to probe new fundamental physics in current and future 
observations, even if the scale of the new fundamental physics is close to the Planck scale. The 
point is that if the universe in fact underwent a period of inflation, then trans-Planckian physics 
leaves an imprint on the spectrum of fluctuations. The exponential expansion of space amplifies 
the wavelength of the perturbations to observable scales. At the present time, it is our ignorance 
about quantum gravity which prevents us from making any specific predictions. For example, we 
do not understand string theory in time-dependent backgrounds sufficiently well to be able to at 
this time make any predictions for observations. 

Singularity Problem 

The next problem is the "singularity problem". This problem, one of the key problems of 
Standard Cosmology, has not been resolved in models of scalar field-driven inflation. 

As follows from the Penrose-Hawking singularity theorems of General Relativity (see e.g. 
[ p9| ] for a textbook discussion), an initial cosmological singularity is unavoidable if space-time is 
described in terms of General Relativity, and if the matter sources obey the weak energy conditions. 
Recently, the singularity theorems have been generalized to apply to Einstein gravity coupled to 
scalar field matter, i.e. to scalar field-driven inflationary cosmology [90]. It is shown that in this 
context, a past singularity at some point in space is unavoidable. 

In the same way that the appearance of an initial singularity in Standard Cosmology told us that 
Standard Cosmology cannot be the correct description of the very early universe, the appearance 
of an initial singularity in current models of inflation tell us that inflationary cosmology cannot 
yield the correct description of the very, very early universe. At sufficiently high densities, a 
new description will take over. In the same way that inflationary cosmology contains late-time 
standard cosmology, it is possible that the new cosmology will contain, at later times, inflationary 
cosmology. However, one should keep an open mind to the possibility that the new cosmology will 
connect to present observations via a route which does not contain inflation, a possibility explored 
later in this section. 
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Breakdown of Validity of Einstein Gravity 
The Achilles heel of scalar field-driven inflationary cosmology is, however, the use of intuition 
from Einstein gravity at energy scales not far removed from the Planck and string scales, scales 
where correction terms to the Einstein-Hilbert term in the gravitational action dominate and where 



intuition based on applying the Einstein equations break down (see also [ ]91| ] for arguments along 
these lines). 

All approaches to quantum gravity predict correction terms in the action which dominate at 
energies close to the Planck scale - in some cases in fact even much lower. Semiclassical gravity 
leads to higher curvature terms, and may (see e.g. [[)2|, ^3|]) lead to bouncing cosmologies without a 
singularity). Loop quantum cosmology leads to similar modifications of early universe cosmology 
(see e.g. [ p4| ] for a recent review). String theory, the theory we will focus on in the following 
sections, has a maximal temperature for a string gas in thermal equilibrium [f>5|], which may lead 



to an almost static phase in the early universe - the Hagedorn phase Q96Q . 

Common to all of these approaches to quantum gravity corrections to early universe cosmology 
is the fact that a transition from a contracting (or quasi-static) early universe phase to the rapidly 
expanding radiation phase of standard cosmology can occur without violating the usual energy 
conditions for matter. In particular, it is possible (as is predicted by the string gas cosmology 
model discussed below) that the universe in an early high temperature phase is almost static. This 
may be a common feature to a large class of models which resolve the cosmological singularity. 

Closely related to the above is the "cosmological constant problem" for inflationary cosmol- 
ogy. We know from observations that the large quantum vacuum energy of field theories does 
not gravitate today. However, to obtain a period of inflation one is using precisely the part of the 
energy-momentum tensor of the inflaton field which looks like the vacuum energy. In the absence 
of a convincing solution of the cosmological constant problem it is unclear whether scalar field- 
driven inflation is robust, i.e. whether the mechanism which renders the quantum vacuum energy 
gravitationally inert today will not also prevent the vacuum energy from gravitating during the 
period of slow-rolling of the inflaton field 10 . 

4.2 String Gas Cosmology 

An immediate problem which arises when trying to connect string theory with cosmology is 
the dimensionality problem. Superstring theory is perturbatively consistent only in ten space-time 
dimensions, but we only see three large spatial dimensions. The original approach to addressing this 
problem was to assume that the six extra dimensions are compactified on a very small space which 
cannot be probed with our available energies. However, from the point of view of cosmology, it is 
quite unsatisfactory not to be able to understand why it is precisely three dimensions which are not 
compactified and why the compact dimensions are stable. Brane world cosmology [^] provides 
another approach to this problem: it assumes that we live on a three-dimensional brane embedded 
in a large nine-dimensional space. Once again, a cosmologically satisfactory theory should explain 
why it is likely that we will end up exactly on a three-dimensional brane (for some interesting work 
addressing this issue see [p9|, [Tooi[ToT| ]). 



I0 Note that the approach to addressing the cosmological constant problem making use of the gravitational back- 
reaction of long range fluctuations (see [ p7| ] for a summary of this approach) does not prevent a long period of inflation 
in the early universe. 
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Finding a natural solution to the dimensionality problem is thus one of the key challenges 
for superstring cosmology. This challenge has various aspects. First, there must be a mechanism 
which singles out three dimensions as the number of spatial dimensions we live in. Second, the 
moduli fields which describe the volume and the shape of the unobserved dimensions must be 
stabilized (any strong time-dependence of these fields would lead to serious phenomenological 
constraints). This is the moduli problem for superstring cosmology. As mentioned above, solving 
the singularity problem is another of the main challenges. These are the three problems which 
string gas cosmology [96, 102, 103| ] explicitly addresses at the present level of development 11 . 



In the absence of a non-perturbative formulation of string theory, the approach to string cos- 
mology which we have suggested, string gas cosmology [96, 102| , |103| ] (see also [105[ for early 



work, and [ |106| , |107 , |108| , |109[ ] for reviews), is to focus on symmetries and degrees of freedom 
which are new to string theory (compared to point particle theories) and which will be part of a 
non-perturbative string theory, and to use them to develop a new cosmology. The symmetry we 
make use of is T-duality, and the new degrees of freedom are string winding modes and string 
oscillatory modes. 

We take all spatial directions to be toroidal, with R denoting the radius of the torus. Strings 
have three types of states: momentum modes which represent the center of mass motion of the 
string, oscillatory modes which represent the fluctuations of the strings, and winding modes count- 
ing the number of times a string wraps the torus. Both oscillatory and winding states are special to 
strings. Point particle theories do not contain these modes. 

The energy of an oscillatory mode is independent of R, momentum mode energies are quan- 
tized in units of l/R, i.e. 

E n = nj, (4.2) 
whereas the winding mode energies are quantized in units of R, i.e. 

E m = mR, (4.3) 

where both n and m are integers. The energy of oscillatory modes does not depend on R. 

The T-duality symmetry is the invariance of the spectrum of string states under the change 

R^ l/R (4.4) 

in the radius of the torus (in units of the string length l s ). Under such a change, the energy spectrum 
of string states is not modified if winding and momentum quantum numbers are interchanged 

(n,m) — ► (m,n) . (4.5) 

The string vertex operators are consistent with this symmetry, and thus T-duality is a symmetry of 
perturbative string theory. Postulating that T-duality extends to non-perturbative string theory leads 
[ |110| ] to the need of adding D-branes to the list of fundamental objects in string theory. With this 
addition, T-duality is expected to be a symmetry of non-perturbative string theory. Specifically, 



1 1 See also [ 104 ] for an alternative approach within string theory to explaining why only three spatial dimensions are 
large. 
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Figure 4: Sketch (based on the analysis of of the evolution of temperature T as a function of the radius 
R of space of a gas of strings in thermal equilibrium. The top curve is characterized by an entropy higher 
than the bottom curve, and leads to a longer region of Hagedorn behaviour. 



T-duality will take a spectrum of stable Type IIA branes and map it into a corresponding spectrum 



of stable Type IIB branes with identical masses [111]. 



Since the number of string oscillatory modes increases exponentially as the string mode energy 
increases, there is a maximal temperature of a gas of strings in thermal equilibrium, the Hagedorn 
temperature Th If we imagine taking a box of strings and compressing it, the temperature 
will never exceed Th- In fact, as the radius R decreases below the string radius, the temperature 
will start to decrease, obeying the duality relation @] 

T(R) = T(l/R). (4.6) 

This argument shows that string theory has the potential of taming singularities in physical observ- 
ables. Figure 4 provides a sketch of how the temperature T changes as a function of R. 

If we imagine that there is a dynamical principle that tells us how R evolves in time, then 
Figure 2 can be interpreted as depicting how the temperature changes as a function of time. If R 
is a monotonic function of time, then two interesting possibilities for cosmology emerge. If In/? 
decreases to zero at some fixed time (which without loss of generality we can call t = 0), and 
continues to decrease, we obtain a temperature profile which is symmetric with respect to t = 
and which (since small R is physically equivalent to large R) represents a bouncing cosmology (see 



[ |1 12| ] for a concrete recent realization of this scenario). If, on the other hand, it takes an inifinite 
amount of time to reach R = 0, an emergent universe scenario [ ]l 13| ] is realized. 

It is important to realize that in both of the cosmological scnearios which, as argued above, 
seem to follow from string theory symmetry considerations alone, a large energy density does not 
lead to rapid expansion in the Hagedorn phase, in spite of the fact that the matter sources we are 
considering (namely a gas of strings) obey all of the usual energy conditions discussed e.g. in 
[|89[]). These considerations are telling us that intuition drawn from Einstein gravity will give us a 
completely incorrect picture of the early universe. 

Any physical theory requires both a specification of the equations of motion and of the initial 
conditions. We assume that the universe starts out small and hot. For simplicity, we take space 
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to be toroidal, with radii in all spatial directions given by the string scale. We assume that the 
initial energy density is very high, with an effective temperature which is close to the Hagedorn 
temperature, the maximal temperature of perturbative string theory. 

In this context, it was argued [ |96| ] that in order for spatial sections to become large, the wind- 
ing modes need to decay. This decay, at least on a background with stable one cycles such as a 
torus, is only possible if two winding modes meet and annihilate. Since string world sheets have 
measure zero probability for intersecting in more than four space-time dimensions, winding modes 
can annihilate only in three spatial dimensions (see, however, the recent caveats to this conclusion 



based on the work of [11141, 115]). Thus, only three spatial dimensions can become large, hence 



explaining the observed dimensionality of space-time. As was shown later [ |103| ], adding branes to 
the system does not change these conclusions since at later times the strings dominate the cosmo- 
logical dynamics. Note that in the three dimensions which are becoming large there is a natural 



mechanism of isotropization as long as some winding modes persist [ 1 16]. 

Some of the above heuristic arguments can be put on a more firm mathematical basis, albeit in 
the context of a toy model, a model consisting of a classical background coupled to a gas of strings. 
From the point of view of rigorous string theory, this separation between classical background and 
stringy matter is not satisfactory when dealing with very early times when the typical length scale 
might be the string scale.. However, in the absence of a non-perturbative formulation of string 
theory, at the present time we are forced to make this separation. Note that this separation between 
classical background geometry and string matter is common to all current approaches to string 
cosmology. 

The background is described by dilaton gravity. The dilaton must be included since it arises 
in string theory at the same level as the graviton, and also (in the context of string gas cosmology) 
because it and not Einstein gravity is consistent with the T-duality symmetry. Note, however, that 
the background dynamics inevitably drives the system into a parameter region where the dilaton is 
strongly coupled and hence beyond the region of validity of the approximations made. 

The action for dilaton gravity coupled to a matter action S m is 

X - 2 J d 10 xV=ge- 2 * [R + 4d^d^] +S m , (4.7) 



2k 

where g is the determinant of the metric, R is the Ricci scalar, (p is the dilaton, and K is the reduced 
gravitational constant in ten dimensions. The metric appearing in the above action is the metric in 
the string frame. 

For a homogeneous and isotropic metric 

ds 2 = dt 2 -a{t) 2 dx 2 , (4.8) 



the resulting equations of motion in the string frame are [102] (see also [1 17]) 

-dX 2 + (p 2 = e^E (4.9) 

X-<pX = (4.10) 

(p-dX 2 = -e^E, (4.11) 
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where E and P denote the total energy and pressure, respectively, d is the number of spatial dimen- 
sions, and we have introduced the logarithm of the scale factor 

X(t) = log(a(0) (4.12) 

and the rescaled dilaton 

(p = 2(j)-dX. (4.13) 

The second of these equations indicates that a gas of strings containing both stable winding 
and momentum modes will lead to the stabilization of the radius of the torus: windings prevent 
expansion, momenta prevent the contraction. The right hand side of the equation can be interpreted 
as resulting from a confining potential for the scale factor. One of the key issues when dealing 
with theories with extra dimensions is the question of how the size and shape moduli of the extra- 
dimensional spaces are stabilized. String gas cosmology provides a simple and string-specific 



mechanism to stabilize most of these moduli. This topic will not be reviewed here (see [ ]108| , |109| ] 
for recent reviews). The outstanding issue is how to stabilize the dilaton. 

Note that the dilaton is evolving at the time when the radius of the torus is at the minimum 
of its potential. For the branch of solutions we are considering, the dilaton is increasing as we 
go into the past. At some point, therefore, it becomes greater than zero. At this point, we enter 



the region of strong coupling. As already discussed in [118], a different dynamical framework 



is required to analyze this phase. In particular, the fundamental strings are no longer the lightest 



degrees of freedom. We will call this phase the "strongly coupled Hagedorn phase" [119] for which 
we lack an analytical description. Since the energy density in this phase is of the string scale, the 
background equations should also be very different from the dilaton gravity equations used above. 
In the following, we assume that the dilaton is frozen in the strongly coupled Hagedorn phase. This 
could be a consequence of S-duality (see e.g. [ ]120| ]). 



4.3 String Gas Cosmology and Structure Formation 

The following are key aspects of the string gas cosmology background which emerge from the 
previous discussion. First, in thermal equilibrium at the string scale (R ~ l s ), the self-dual radius, 
the number of winding and momentum modes are equal. Since winding and momentum modes 
give an opposite contribution to the pressure, the pressure of the string gas in thermal equilibrium 
at the self-dual radius will vanish. From the dilaton gravity equations of motion ( |4.9| - [4.1 1[ ) it then 
follows that a static phase A = will be a fixed point of the dynamical system. This phase is the 
Hagedorn phase. 

On the other hand, for large values of R in thermal equilibrium the energy will be exclusively 
in momentum modes. These act as usual radiation. Inserting the radiative equation of state into the 



above equations ( |4.9| - 14.11| ) it follows that the source in the dilaton equation of motion vanishes 
and the dilaton approaches a constant as a consequence of the Hubble damping term in its equation 
of motion. Consequently, the scale factor expands as in the usual radiation-dominated universe. 
The transition between the Hagedorn phase and the radiation-dominated phase with fixed dilaton 



is achieved via the annihilation of winding modes, as studied in detail in [121]. The main point 
is that, starting in a Hagedorn phase, there will be a smooth transition to the radiation-dominated 
phase of standard cosmology with fixed dilaton. 
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Our new cosmological background is obtained by following our currently observed universe 
into the past according to the string gas cosmology equations. The radiation phase of standard 
cosmology is unchanged. In particular, the dilaton is fixed in this phase 12 . However, as the tem- 
perature of the radiation bath approaches the Hagedorn temperature, the equation of state of string 
gas matter changes. The equation of state parameter w = P/E decreases towards a pressureless 
state and the string frame metric becomes static. Note that, in order for the present size of the 
universe to be larger than our current Hubble radius, the size of the spatial sections in the Hagedorn 
phase must be at least 1mm 13 . We have entered the Hagedorn phase. 

As we go back in time in the Hagedorn phase, the dilaton increases. At the time t c when 
the dilaton equals zero, a second transition occurs, the transition to a "strongly coupled Hagedorn 



phase" (using the terminology introduced in [1 19]). We take the dilaton to be fixed in this phase. In 



this case, the strongly coupled Hagedorn phase may have a duration which is very long compared 
to the Einstein frame Hubble time immediately following tR. It is in this cosmological background 
that we will study the generation of fluctuations. We will denote the time when the Einstein frame 
Hubble radius reaches a minimum by tR, to evoque the analogy with the time of reheating in infla- 
tionary cosmology. The end of the Hagedorn phase (the time when the radiation phase of standard 
cosmology begins) is slightly later. 

It is instructive to compare the background evolution of string gas cosmology with the back- 
ground of inflationary cosmology. Figure 5 is a sketch of the space-time evolution in string gas 
cosmology. For times t < tR, we are in the static Hagedorn phase and the Hubble radius is infinite. 
For t > tR, the Einstein frame Hubble radius is expanding as in standard cosmology. To understand 
why string gas cosmology can lead to a causal mechanism of structure formation, we must com- 
pare the evolution of the physical wavelength corresponding to a fixed comoving scale with that of 
the Einstein frame Hubble radius H~ l {t). Recall that the Einstein frame Hubble radius separates 
scales on which fluctuations oscillate (wavelengths smaller than the Hubble radius) from wave- 
lengths on which the fluctuations are frozen in and cannot be effected by microphysics. Causal 
microphysical processes can generate fluctuations only on sub-Hubble scales 14 . The key point is 
that for t < tj{k), the fluctuation mode k is inside the Hubble radius, and thus a causal generation 
mechanism for fluctuations is possible. 

In contrast, in inflationary cosmology (Figure 1) the Hubble radius is constant during inflation 
(t < tR, where here tR is the time of inflationary reheating), whereas the physical wavelength corre- 
sponding to a fixed comoving scale expands exponentially. Thus, as long as the period of inflation 
is sufficiently long, all scales of interest for current cosmological observations are sub-Hubble at 
the beginning of inflation. 

There are both important similarities and key differences between the structure formation 
mechanisms in inflationary cosmology and string gas cosmology. In both cases, scales are sub- 
Hubble during the early stages, thus allowing for a causal generation mechanism. Also, in both 



12 The dilaton comes to rest, but it is not pinned to a particular value by a potential. Thus, in order to obtain consis- 
tency with late time cosmology, an additional mechanism operative at late times which fixes the dilaton is required. 



13 



How to obtain this initial size starting from string-scale initial conditions constitutes the entropy problem of our 



scenario. A possible solution making use of an initial phase of bulk dynamics is given in [122]. 

14 Matter which couples minimally to gravity in the string frame is prevented from oscillating on scales larger than 
the Einstein frame Hubble radius by the dilaton friction term in its equation of motion. 
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Figure 5: Space-time diagram (sketch) showing the evolution of fixed comoving scales in string gas cos- 
mology. The vertical axis is time, the horizontal axis is physical distance. The solid curve represents the 
Einstein frame Hubble radius H~ l which shrinks abruptly to a microphysical scale at tR and then increases 
linearly in time for t > tR. Fixed comoving scales (the dotted lines labeled by k\ and h£) which are currently 
probed in cosmological observations have wavelengths which are smaller than the Hubble radius before tR. 
They exit the Hubble radius at times ?, (£) just prior to tR, and propagate with a wavelength larger than the 
Hubble radius until they reenter the Hubble radius at times tf(k). 

cases the fluctuations evolve on super-Hubble scales for a long time in the radiation phase of stan- 
dard cosmology, thus leading to their squeezing, which in turn leads to the phase coherence of the 
fluctuations which generate the acoustic oscillations in the angular power spectrum of the CMB. 
However, the actual generation mechanism for fluctuations is completely different. In inflation- 
ary cosmology, any thermal fluctuations present before the onset of inflation are red-shifted away, 
leaving us with a quantum vacuum state, whereas in the quasi-static Hagedorn phase of string gas 
cosmology matter is in a thermal state. Hence, whereas in inflationary cosmology the fluctua- 
tions originate as quantum vacuum perturbations, in string gas cosmology the inhomogeneities are 
created by the thermal fluctuations of the string gas. 



As we have shown in [ |123| , |124] , |125| ], string thermodynamical fluctuations in the Hagedorn 
phase of string gas cosmology yield an almost scale-invariant spectrum of both scalar and tensor 
modes. This result stems from the holographic scaling of the specific heat Cy(R) (evaluated for 
fixed volume) as a function of the radius R of the box 

C V {R)~R 2 . (4.14) 



As derived in [126], this result holds true for a gas of closed strings in a space-time in which the 
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three large spatial dimensions are compact (see []129|, |130|] for recent papers emphasizing the role 



of holography). The scaling ( |4. 14{ ) is an intrinsically stringy result: thermal fluctuations of a gas of 
particles would lead to a very different scaling. 

Since the primordial perturbations in our scenario are of thermal origin (and there are no non- 
vanishing chemical potentials), they will be adiabatic. The spectrum of scalar metric fluctuations 



has a slight red tilt. As a distinctive feature [127], our scenario predicts a slight blue tilt for the 



spectrum of gravitational waves. The red tilt for the scalar modes is due to the fact that the temper- 
ature when short wavelength modes exit the Hubble radius is slightly lower than the temperature 
when longer wavelength modes exit. The gravitational wave amplitude, in contrast, is determined 
by the pressure. Since the pressure is closer to zero the deeper in the Hagedorn phase we are, a 
slight blue tilt for the tensor fluctuations results. These results are explained in more detail in a 



recent review []109p an d i n tne original references. Here, we will only focus on some of the key 
steps. 

Our approximation scheme for computing the cosmological perturbations and gravitational 
wave spectra from string gas cosmology is as follows (the analysis is similar to how the calculations 



were performed in Q54| , |55[ ] in the case of inflationary cosmology). For a fixed comoving scale k 
we follow the matter fluctuations until the time tj{k) shortly before the end of the Hagedorn phase 
when the scale exits the Hubble radius 15 At that time, we use the Einstein constraint equations 
(discussed below) to compute the values of <&(k) and h{k) (h is the amplitude of the gravitational 
wave tensor), and then we propagate the metric fluctuations according to the standard gravitational 
perturbation equations until scales re-enter the Hubble radius at late times. 

The first key point is to show how the scalar and tensor metric fluctuations can be extracted 
from knowledge of the energy-momentum tensor of the string gas. Working in conformal time v\ 
and in the longitudinal gauge for the scalar metric fluctuations, the metric of a homogeneous and 
isotropic background space-time perturbed by linear cosmological perturbations and gravitational 
waves can be written in the form 

ds 2 = a 2 {ri) { (1 + 2<!>)dri 2 - [(1 - 2*)fy + h^dx'dx-l } . (4.15) 

Here, <I> (which is a function of space and time) describes the scalar metric fluctuations 16 . The 
tensor hij is transverse and traceless and contains the two polarization states of the gravitational 
waves. We have assumed that there is no anisotropic stress. 

Inserting the metric ( |4. 15 ) into the Einstein equations, subtracting the background terms and 



truncating the perturbative expansion at linear order leads to the following system of equations 



- 3^ ( Jt?& + <£') +V 2 <£ = 47iGa 2 8T° 



o 



+ ^ = 4nGa 2 8T°i 



I5 Recall that on sub-Hubble scales, the dynamics of matter is the dominant factor in the evolution of the system, 
whereas on super-Hubble scales, matter fluctuations freeze out and gravity dominates. Thus, it is precisely at the time 
of Hubble radius crossing that we must extract the metric fluctuations from the matter perturbations. Since the concept 
of an energy density fluctuation is gauge-dependent on super-Hubble scales, one cannot extrapolate the matter spectra 
to larger scales as was done in Section 3 of [ |l28| ]. 

I6 To avoid confusion with the dilaton <j>, we in this section denote the relativistic gravitational potential by <!>. 
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2Jif + Jf 2 <J> + 3^<£ + <£ 



J^' + -Jif 2 
2 



-4KGa 2 8rj, 
for i^j. 



(4.16) 



In the Hagedorn phase, these equations simplify substantially and allow us to extract the scalar 
and tensor metric fluctuations individually. Replacing comoving by physical coordinates, we obtain 
from the 00 equation 



V 2 <£ = 4jiGST° 



and from the i ^ j equation 



W 2 h, 



-AnGST*} . 



(4.17) 



(4.18) 



The above equations ( 117 ) and (118) allow us to compute the power spectra of scalar and 
tensor metric fluctuations in terms of correlation functions of the string energy-momentum tensor. 
Since the metric perturbations are small in amplitude we can consistently work in Fourier space. 
Specifically, 

{\<S>{k)\ 2 ) = 16k 2 G 2 k- 4 (8T° (k)8T° Q (k)), (4.19) 
where the pointed brackets indicate expectation values, and 



(\h(k)\ 2 ) = \6x 2 G 2 k- 4 (ST'j(k)ST'j(k)), 



(4.20) 



where on the right hand side of ( |120[ ) we mean the average over the correlation functions with 

i + h 

The second key step is to compute the matter fluctuations in the Hagedorn phase. Since this 
phase is dominated by the gas of strings, fluctuations in our scenario are the thermal fluctuations of 
a string gas. We will consider a gas of closed strings in a compact space, i.e. our three-dimensional 
space is considered to be large but compact. Specifically, it is important to have winding modes in 
the spectrum of string states. 

General thermodynamical relations allow us to compute the matter correlation functions in 
terms of the the specific heat Cy and the pressure p. The result for the energy density fluctuation is 



(4.21) 



where F is the free energy of the string gas and j3 is the inverse temperature. The off-diagonal 
pressure fluctuations, in turn, are given by 



i \2 



(S7V) = (T'j ) - (T'j) 



1 dF 



f$R 3 d InR \ R 3 d InR J fiR 2 dR 



1 dp 



(4.22) 



Now, we apply these relations to the thermodynamics of strings. In [ ]126| ], the thermodynami- 
cal properties of a gas of closed strings in a toroidal space of radius R were computed. To compute 
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the fluctuations in a region of radius R which forms part of our three-dimensional compact space, 
we will apply the results of [126] to a box of strings in a volume V = R . 



The starting point of the computation is the formula for the density of states Q.(E,R) which 
then determines the cpecific heat and the pressure. The result for the specific heat is 



C v 



R 2 /£ 3 



T(\-T/T H ) 



(4.23) 



The 'holographic' scaling Cy (R) ~ R 2 is responsible for the overall scale-invariance of the spectrum 
of cosmological perturbations. The factor (1 — T /Tjj) in the denominator is responsible for giving 
the spectrum a slight red tilt. 
For the pressure, we obtain 



p{E,R) ~ nfjTfj 



2 {\-T/T H ) 
3 @R 



In 



R 2 (l-T/T H ) 



(4.24) 



which immediately yields 



<^>.I™ln 2 



(4.25) 



Note that the factor (1 — T /Th) (which appears in the numerator of the key term in the expression 
for the microcanonical partition function) is in the numerator, whereas it was in the denominator 
in the expression for the specific heat. The reason is that in deriving the specific heat from the 
microcanonical partition function, a temperature derivative was taken, but not so in deriving the 
pressure. This leads to the slight blue tilt of the spectrum of gravitational waves, characteristic of 
our proposed structure formation scenario. As mentioned earlier, the physical reason for this blue 
tilt is that larger wavelength modes exit the Hubble radius deeper in the Hagedorn phase where the 
pressure is smaller and thus the strength of the tensor modes is less. 

The third key step in the analysis is to compute the power spectra for the scalar and tensor 
modes based on the earlier results. The power spectrum of scalar metric fluctuations is given by 



(4.26) 



SG 2 k~ l < \8p(k)\ 2 > 



%G l k l < (8M) Z > R 
SG 2 k- 4 <(Sp) 2 > R , 



where in the first step we have used ( 4.19| ) to replace the expectation value of | <!>(&) | 2 in terms of 
the correlation function of the energy density, and in the second step we have made the transition 



to position space (note that k = R *). 



According to (4.21), the density correlation function is given by the specific heat via T 2 R- 6 C V . 
Inserting the expression from ( |4.23 ) for the specific heat of a string gas on a scale R yields to the 
final result 



Bs,(k) = 8G 2 



4 3 1 - t/t h 



(4.27) 
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for the power spectrum of cosmological fluctuations. In the above equation, the temperature T is 
to be evaluated at the time tj{k) when the mode k exits the Hubble radius. Since modes with larger 
values of k exit the Hubble radius slightly later when the temperature is slightly lower, a small red 
tilt of the spectrum is induced. The amplitude sis of the power spectrum is given by 

1 



sis 



T/T H 



(4.28) 



Taking the last factor to be of order unity, we find that a string length three orders of magnitude 
larger than the Planck length, a string length which was assumed in early studies of string theory, 
gives the correct amplitude of the spectrum. Thus, it appears that the string gas cosmology structure 
formation mechanism does not have a serious amplitude problem. 

Similarly, we can compute the power spectrum of the gravitational waves and obtain 



P h (k) ~ SG 2 p(l-T/T H )ln 2 



1 



e 2 k 2 



(l-T/Ta 



,-1 



(4.29) 



This shows that the spectrum of tensor modes is - to a first approximation, namely neglecting the 
logarithmic factor and neglecting the k-dependence of T (/,•(&)) - scale-invariant. The k-dependence 
of the temperature at Hubble radius crossing induces a small blue tilt for the spectrum of gravita- 
tional waves. 



Comparing (4.27) and (4.29) we see that the tensor to scalar ratio is suppressed by the factor 
(1 — T /Th) 2 . Given a good understanding of the exit from the Hagedorn phase we would be able 
to compute this ratio as well as the magnitude of the spectral tilts for both scalar and tensor modes. 

4.4 Discussion 

In order to put our structure formation scenario on a firm basis, we need a consistent descrip- 
tion of the Hagedorn phase. The dilaton gravity background allows us to understand the onset of 
the Hagedorn phase (going backwards in time), but since the dilaton blows up, we rapidly leave the 
domain of applicability of the model. 



A background in which our string gas structure formation scenario can be implemented [1 12] 



is the ghost-free and asymptotically free higher derivative gravity model proposed in [ |93| ] given by 
the gravitational action 

S = f d A x^f^F(R) (4.30) 



with 



FIR) = R+Y -%-R&r 



v 2 Yr, 



(4.31) 



where M„ is the string mass scale (more generally, it is the scale where non-perturbative effects 
start to dominate), and the c n are coefficients of order unity. 



As shown in [ ]93[ ] and [ |112| ], this action has bouncing cosmological solutions. If the tem- 
perature during the bounce phase is sufficiently high, then a gas of strings will be excited in this 
phase. In the absence of initial cosmological perturbations in the contracting phase, our string gas 
structure formation scenario is realized. The string network will contain winding modes in the 
same way that a string network formed during a cosmological phase transition will contain infi- 
nite strings. The dilaton is fixed in this scenario, thus putting the calculation of the cosmological 



38 



Topics in Cosmology 



Robert Brandenberger 



perturbations on a firm basis. There are no additional dynamical degrees of freedom compared to 
those in Einstein gravity. The higher derivative corrections to the equations of motion (in partic- 
ular to the Poisson equation) are suppressed by factors of (k/M s ) 2 . Thus, all of the conditions on 
a cosmological backgroound to successfully realize the string gas cosmology structure formation 
scenario are realized. 

5. Conclusions 

These lectures have focused on three topics in theoretical cosmology. The first is the infla- 
tionary universe scenario, the current paradigm of early universe cosmology. Inflation has been 
an extremely successful scenario. It explains why the universe is spatially flat and isotropic to the 
extent it is observed to be, and why it is so large and contains such a large amount of entropy, thus 
resolving several mysteries which the previous paradigm of cosmology, the SBB model, was not 
able to explain. Possibly more importantly, inflation provides a causal mechanism for the origin 
of the fluctuations which are now mapped out to great accuracy by recent cosmological observa- 
tions. Inflation correctly (and fifteen years before the precision observations) predicted the power 
spectrum of CMB anisotropics. 

The theory of cosmological perturbations is the key tool which allows us to take theories of 
the very early universe and calculate predictions for late time cosmology. This theory is applicable 
whatever the paradigm of early universe cosmology might be, and it is the second topic discussed 
in these lectures. 

In spite of the successes of the inflationary universe scenario, key conceptual questions remain, 
making it clear that new input from fundamental physics is required in order to develop a better 
theory of the very early universe. The third topic discussed in these lectures, string gas cosmology, 
is an attempt to explore possible consequences for cosmology of superstring theory, the likely 
candidate for the new physics required for early universe cosmology. In particular, an alternative 
structure formation scenario not requiring inflationary dynamics emerges. 
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